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Abstract 



5-1 

Qh! We consider a dynamical system with state space M, a smooth, compact subset of 



some M. n , and evolution given by T t , xt = T t x, x G M; T t is invertible and the time 
t may be discrete, t G Z, T t = T*, or continuous, i 6 1. Here we show that starting 



with a continuous positive initial probability density p(x, 0) > 0, with respect to dx, the 
3 ! smooth volume measure induced on M by Lebesgue measure on M n , the expectation value 
of log/9(x,t), with respect to any stationary (i.e. time invariant) measure v(dx), is linear 

00 ' 

^| absolutely continuous wrt dx. 

a 



in t, u(logp(x,t)) = v(logp(x, 0)) + Kt. K depends only on v and vanishes when v is 



^ ■ The time evolution of probability measures on the phase space M of a dynamical 

system depends both on the character of the dynamics, assumed here to be given by a one 
parameter group of invertible measurable transformations Tt, and the nature of the initial 
measure. Given a probability measure po on M at time 0, the evolved measure at time t, 
fit, is such that the expectation value of functions (f>(x) is given by 

Pt{4>) = / 4>(x)p t (dx) = / <p(T t x)fj,o(dx), (1) 
J m J 



or in terms of measurable sets A G M 



IH{A) = no(T- t A) 
1 



(2) 



where T_ t A is the set of points y such that T t y belongs to A. 

There will typically be many stationary measures v, ut(dx) = v(dx), for the dynamical 
system. Some are concentrated on the set of fixed points or periodic orbits. There 
can also be stationary measures concentrated on fractal sets of zero Lebesgue measure. 
This happens in particular for generic Anosov systems; for other examples, see [1,2]. In 
addition there may exist a stationary measure absolutely continuous with respect to dx, 
i.e. v(dx) = p(x)dx. In the most familiar examples of the latter situation T t preserves dx, 
as with Hamiltonian flows on energy surfaces or the baker's transformation on the unit 
square, in which case p(x) is constant, i.e. independent of x. 

When T t is sufficiently "chaotic", the fractal and absolutely continuous types of 
measures will generally have very good ergodic properties, mixing or better. For the 
case of an absolutely continuous v, mixing implies that if fio(dx) = p(x,0)dx then 

IhU) = J f(T t x)p(x,0)dx = J f(x)p(x,t)dx J f(x)p(x)dx (3) 

for any bounded measurable f(x). The weak convergence of p(x,t) to p, expressed by 
(3), is clearly compatible with the fact that when T t preserves dx, the Gibbs entropy 
= — f p(x, t) logp(x, t)dx, and indeed any F = f F(p(x, t))dx, is constant in time. 
Unfortunately, it is sometimes thought that this constancy of for Hamiltonian 
evolutions is a manifestation of the conflict between microscopic reversibility and the 
second law of thermodynamics, and that the resolution of this conflict requires at least 
an acceptance of weak convergence as the mathematical expression of the approach to 
equilibrium characteristic of macroscopic irreversibility, and perhaps even necessitates 
changes in the microscopic physical laws, c.f. [3b]. This concern and its proposed 
resolution are based on a misunderstanding of the origin of the observed time asymmetry 
of macroscopic physical systems, which really concerns not probability densities but the 
behavior of individual systems whose microstates x t = T t x are points in a very high 
dimensional phase space M. In fact, the second law refers not to but to an entropy 
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defined for individual macroscopic systems, whose observed irreversible behavior is due 
first and foremost to the large discrepancy between the scale of macroscopic observables 
(which behave irreversibly) and microscopic scales and to the nature of "typical" initial 
conditions for the microstate x of the system, c.f. [3]. 

These conceptual issues are, however, not the main concern of this brief note, even 
though it was motivated by the paper of R. Fox in this issue [4] in which such problems 
are discussed for the baker's transformation. In that paper Fox notes the constancy of 
F when F = logp for this transformation. Here we are concerned with what happens to 
functions of p(x,t) when T t does not preserve dx and v may not be absolutely continuous. 

Let v be a stationary probability measure, and let p^ and p^ be two measures on 
M, evolving according to the dynamics. If p^ is absolutely continuous wrt to p±\ i.e. 
p[ 1 \dx) = g(x,t)p[ 2 \dx), then it follows directly from (2) that 

g(x,t) = g(T- t x,0). (4) 

Suppose that g(x,0) is continuous in x. Then, given a function of g, f(g), integrable 
wrt v, we have v(f(g(x,t))) = v(f(g(x,0))) for all t. Assume now that p^ and p\ 2 ^ are 
themselves absolutely continuous wrt dx, with continuous positive densities pi(x,t) and 
P2(x,t). Then g(x, t) = pi(x, t)/p2(x, t) and 

Hf(g))= [ f{pi(x,t)/p 2 (x,t))u{dx)= Const. (5) 

J M 

Setting f(g) = logg yields 

i/(logpi(x,t)) - v(\ogp 2 (x,t)) = C (6) 

independent of t. Put now p2(x,t) = p(x,t) and pi(x,t) = p(x,t + r). Eq. (5) then 
becomes for all r 

v(\ogp(x,t + T)) - is(\ogp(x,t)) = K(t). (7) 
3 



Noting that K{t\ + r 2 ) = K(t\) + K(t 2 ) we obtain a rather surprising result 

i/(logp(x,r)) = i/(logp(x,0)) + iirT (8) 

with K independent of r. In other words, the average of the log of the density wrt the 
stationary measure v is linear in the time. On the other hand it follows from (6) that the 
growth rate of v(\ogp(x, £)) does not depend on p. Hence K depends only on the dynamics 
T t and the stationary probability measure v. Consequently, we can compute K by taking 
for our initial (unnormalized) density p(x, 0) = 1. We then get 

*-(fu), (9) 

where J(x,i) is the Jacobian of the transformation T_ t , for continuous time and 

K = v{\ogJ{x)) (10) 

where J(x) = J(x, 1), for discrete time. If v is absolutely continuous wrt dx, i.e. u(dx) = 
p(x)dx, then putting p2(x,t) = p(x) and p\ = p in (5) we see that J M [log p{x, t)]p(x)dx is 
independent of t, i.e., K vanishes for such a v. 

In the case of a continuous time evolution given by a (smooth) vector field, x = v(x), 
the right side of (9) is just v{— V • v). Eqs. (8) and (9) can then also be obtained directly 
for a smooth, positive p{x, 0) by starting with the continuity equation 

dp(x,t) 



dt 

We then find 



-V-(pv(aO). (11) 



K = ^- I logp(x, t)u(dx) = — [ p X V ■ (pv)v(dx) 
dt J M J M 

= -/ [V- v+ (Vlogp) • v]u(dx) 

JM 

On the other hand, the time derivative of pt{4>) is? f° r an Y smooth 4>(x), given by 



(12) 



1^(0) = (v- V0). (13) 



Hence, by the stationarity of v, v(v ■ V</>) = and so the second term in the square bracket 
in (12) vanishes, yielding explicitly 

K = — z/(V ■ v). (14) 

Eqs. (7) and (13) are to be compared with what happens to the rate of change of the 
Gibbs entropy S^, for p t (dx) = p(x,t)dx. A straightforward computation gives 

d d f f 

jj.Sn = -— J p(x, t) log p(x,t)dx = J (V-v)p(x,t)dx = /x t (V-v). (15) 

Sfj, has been of much interest recently in connection with "thermostatted" nonequilibrium 
systems [1,2,5]. Under suitable conditions on T t , it can be shown that p,t(dx) — > v±(dx) 
with v+ an SRB measure [1,2]. In such cases 

4 t t -^±(V-v) (16) 

at t^±oo 

with z/ + (V ■ v) < 0. The equality holds if and only if v + is absolutely continuous wrt dx, 
i.e. v+ = p + (x)dx. On the other hand when T t is "time reversible" in the sense that there 
exists a transformation R on M, preserving dx, such that R 2 = I and RT t x = T_ t Rx, 
then [1,2] 

K+ = -z/+(V- v) = i/_(V- v) = -K-. (17) 
Thus, writing S±(t) = —v±(log p(x,t)) we have in this case that 

S li (t)~S±(t) for t^±oo 

with 

S±(t) = S±(Q)TK+t. 

As an illustrative example consider a flow on a circle, with v(x) = —sin a; + ui. 
Here x G [— tt, tt] with periodic boundary conditions and a; is a constant. This example 
corresponds to a particle moving in the plane with velocity u under the action of an electric 
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field E and a magnetic field h perpendicular to the plane. The speed |u| is kept equal to one 
by a Gaussian thermostat [2,5]; x is the angle between the velocity u and E and uj ~ h/\E\. 
(This flow is time reversible, with R given by reflection through tt/2, the minimum of v.) 
For \u\ < 1, u±(dx) are delta functions at x± = arcsincu with \x+\ < \x-\. We clearly 
have K + = cosx + = y/1 — uj 2 = —K_ > 0. For |u;| > 1 there is a unique stationary state, 
v{dx) = p(x)dx, with p(x) proportional to l/\v (x)\ and K = on general grounds as well 
as by explicit computation. At \u\ = 1, x + = X-, u + = v- = v with K = so K is 
continuous in uj. 

Another observation which follows from (5) is that for an absolutely continuous v, 
with density p(x), 

B p = f -l\Pp( x )dx (18) 

Jm P\ x ) 

is independent of t. For p = 1 (17) is just the L\ distance between p t and v\ since 
M is compact, f M dx = \M\ < oo, (17) also implies, by the Schwartz inequality, that 
Im \p( x 't) ~ p{x)\ 2 dx > \M\~ 1 Bf > unless p = p, and a similar statement is true of the 
higher norms. Thus there can be no convergence to zero of the L2 and higher norms of 
p(x,t) - p(x). 

We conclude by noting that the long time behavior of p(x,t) was discussed in [6] for 
hyperbolic maps. It was explained there that conditional probability densities induced 
by on the unstable manifolds converge, as t — > 00, pointwise with their derivatives to 
the corresponding densities given by v. Along stable directions, however, the densities 
p(x,t) are extremely irregular, as might be suggested by the preservation of the integrals 
discussed above. 
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